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1 Introduction 

The relation between the two most important in mathematical physics first order 
systems of partial differential equations is among those topics which attract 
attention because of their general, even philosophical significance but at the 
same time do not offer much for the solution of particular problems concerning 
physical models. The Maxwell equations can be represented in a Dirac like form 
in different ways (e.g., ||, ||). Solutions of Maxwell's system can be related 
to solutions of the Dirac equation through some nonlinear equations (e.g., ]Tl[|). 
Nevertheless, in spite of these significant efforts there remain some important 
conceptual questions. For example, what is the meaning of this close relation 
between the Maxwell system and the Dirac equation and how this relation is 
connected with the wave-particle dualism. In the present article we propose a 
simple equality involving the Dirac operator and the Maxwell operators (in the 
sense which is explained below). This equality establishes a direct connection 
between solutions of the two systems and moreover, we show that it is valid 
when a quite natural relation between the frequency of the electromagnetic wave 
and the energy of the Dirac particle is fulfilled. Our analysis is based on the 
quaternionic form of the Dirac equation obtained in j^j] and on the quaternionic 
form of the Maxwell equations proposed in (see also ||). In both cases the 
quaternionic reformulations are completely equivalent to the traditional form of 
the Dirac and Maxwell systems. 
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2 Preliminaries 



The algebra of complex quaternions is denoted by H(C). Each complex quater- 
nion a is of the form a = X)fe=o °fc*fc where {ak} C C, io is the unit and 
{ik\ k = 1, 2, 3} are the quaternionic imaginary units: 



' A-' 



ioh = hio = ik, k = 1,2,3; 



lll 2 = -llH = «3! «2«3 



»3«1 



"ll«3 = «2- 



The complex imaginary unit i commutes with i k , k = 0,3. 

We will use the vector representation of complex quaternions: a = Sc(a) + 
Vec(a) , where Sc(a) = eto and Vec(et) = ~a*=X)fe=i a kik- That is each complex 
quaternion is a sum of its scalar part and its vector part. Complex vectors we 
identify with complex quaternions whose scalar part is equal to zero. In vector 
terms, the multiplication of two arbitrary complex quaternions a and b can be 
written as follows: 



a ■ b = ao^o— < ~ct, b > + if x b + clq b + bo a 



where 



<~ct,b >:= ^^atbk £ 



and 



x b 1 := 



k=l 



b\ b 2 



«3 

«3 
&3 



We shall consider continuously differentiable H(C)— valued functions depend- 
ing on three real variables x = (%i,X2,Xs). On this set the well known (see, 
e.g., Ip, B, 0) Moisil-Theodoresco operator is defined by the expression 



d 



D := y^ikdk, where d k = ^ — . 
^ ox k 



k=l 



The action of the operator D on an H(C)— valued function / can be written in 
a vector form: 
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D/ = -div / +grad/ + rot /. (1) 

That is, Sc(Df) = — div / and Vec(Df) = grad /q + rot / . In a good number 
of physical applications (see M and ||) the operators D a — D + M a and 
D- a = D — M a are needed, where a is a complex quaternion and M a denotes 
the operator of multiplication by a from the right-hand side: M a f = f -a. Here 
we will be interested in two special cases when a is a scalar, that is a = ao or 
when a is a vector a = ~ct. The first case corresponds to the Maxwell equations 
and the second to the Dirac equation. 

3 The Dirac equation 

Consider the Dirac equation in its covariant form 

3 

& + Y] Ikdk) + imc)$(t, x) = 0. 
c f— ' 

fe=i 

For a wave function with a given energy we have x) = q(x)e l ^ t , where q 
satisfies the equation 

,i£ „ imc, , , 

(^70 + 52^ + ^-)^)=°- ( 2 ) 
fe=i 

Denote 

. c 3 

Let us introduce an auxiliary notation / := /(t, a:i, X2, — xq). The transfor- 
mation which allows us to rewrite the Dirac equation in a quaternionic form we 
denote as A and define in the following way Qj. A function $ : R 3 — > C 4 is 
transformed into a function F : R 3 — > H(C) by the rule 

F = A[9] = -(-($1 - *2)«o + i(*o - $3)*i " 

($0 + $3)i2 + i(*l + $2)«3)- 

The inverse transformation A^ 1 is defined as follows 

® = A- 1 {F] = 



(-iFi - F 2 , -F - iF 3> F - iF 3 , iJi - F 2 ). 
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Let us present the introduced transformations in a more explicit matrix form 
which relates the components of a C 4 -valued function <£> with the components 
of an H(C)-valued function F: 



F = A[$] = 





f ° 


-1 


1 









( $0 


\ 


1 


i 








—i 










2 


-1 








-1 














i 


i 





I 









and 



/ 





—i 


-1 





\ 


/ 


% 




-1 








—i 






h 






1 








—i 










\ 





i 


-1 





/ 


V 


F 3 


) 



$ = A- 1 \F] = 



We have the following important equality 



(3) 



where ~ct 



(||) we can consider the equation 



(i—ii + raci-i). This equality shows us that instead of equation 



(4) 



and the relation between solutions of (||) and (|j) is established by means of the 
invertible transformation A: f = Aq. 



4 The Maxwell equations 

We will consider the time-harmonic Maxwell equations for a sourceless isotropic 
homogeneous medium 



rot ~H = —iujel^, 
rot ~~E — iujfiTf , 
div ~H = 0, 
div it = 0. 



(5) 
(6) 
(7) 
(8) 



Here lo is the frequency, e and /i are the absolute permittivity and permeability 
respectively, e = £q£ t and /i = /ioAV, where Sq and /iq arc the corresponding 
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parameters of a vacuum and s r , [i r are the relative permittivity and permeability 
of a medium. 

Taking into account ([!]) we can rewrite this system as follows 

Dif = iuj0, (9) 

DT? = -iwet. (10) 

This pair of equations can be diagonalizcd in the following way || (see also ||). 
Denote 

1$ := -iive + nlf (11) 

and 

~$ := iujet + nil ", (12) 

where n := to^/ejl = ^^/£rfJ-r is the wave number. Applying the operator D to 
the functions ~Tf? and ip one can see that satisfies the equation 

(D-K)-ft = 0, (13) 

and ip satisfies the equation 

(D + k)^ = 0. (14) 
Solutions of ( |l3| ) and (|l4|) are called the Beltrami fields (see, e.g., JnJ). 

5 The relation 

In the preceding sections it was shown that the Dirac equation (^) is equivalent 
to the equation D-&f = with ~ct = — 4(i-«i+mc«2) and the Maxwell equations 
®"(^|) are equivalent to the pair of quaternionic equations D_ K Tff = and 
D K -0=0. Now we will show a simple relation between these objects. Suppose 
that 

k 2 = ~ct 2 . (15) 
Let us introduce the following operators of multiplication 



2k 

It is easy to verify that they are mutually complementary and orthogonal pro- 
jection operators, and the following equality is valid || 

D-g = P+D K + P-D_ K . (16) 
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Moreover, as P ± commute with D± K , we obtain that any solution of (Q) is 
uniquely represented as follows 

where ip and ip are solutions of (|l^) and (|l4|) respectively but in general can be 
full quaternions not necessarily purely vectorial. In particular, we have that 

/ = P + (iiuel! + nit) + p-(-iue~3 + nit) = 
iuje(P + - P-0 + k(P+ + P-0 = 

iuje—> —> 
E ■ o> + kH 

K 

is a solution of (|) if ~E! and H are solutions of 

It should be noticed that ( |l6| ) works in both directions. We have 

D K = P+D^ + P-D--& 

and 

£>_ K = P-D-* + P + D^. 

The fact that the Maxwell system reduces to equations ([l3|) and (|l4]), where the 
functions ~T$ and ip are purely vectorial provokes the natural question whether 
it had any sense to consider full quaternions ip and ip and hence four-component 
vectors E and H or the nature definitely eliminated their scalar parts. Some 
arguments supporting the idea of nonzero scalar parts can be found, for example, 



As we have seen equality ( |16| ) is valid under the condition (|15|). Let us 
analyze this condition. Note that 

a = — < a , a >= — m c ). 



Thus (1l5h has the form 



or equivalently 



(17) 



(Huj) 2 e r (i r = £ 2 — m 2 c 4 
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From this equation in the case e r = [i r = 1, that is for a vacuum, using the well 
known in quantum mechanics relation between the frequency and the impulse: 
hcu = pc we obtain the equality 

£ 2 =p 2 c 2 +m 2 c\ (18) 

In general, if in ( JT^ ) we formally use the de Broglie equality p = Hk, we again 
obtain the fundamental relation (|l8|). 

Thus relation (|l6|) between the Dirac operator and the Maxwell operators 
is valid if the condition (|l7]) is fulfilled which quite surprisingly is in agreement 
with @. 
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